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Field Of The Invention 

The present invention relates generally to a method and apparatus for 
measuring delay in telecommunications signals and is particularly suited to 
determining echo path delay for use with echo cancellers. 

Background 

Echo cancellers for use in telecommunication systems operate by 
developing a model for the echo path impulse response, and then convolving the 
modeled response with an incoming signal in order to produce an estimated echo 
signal. The estimated echo signal is then subtracted from an outbound signal from 
the echo path in order to eliminate echo. A critical parameter required to be 
determined is the time delay of the echo path. For rapid convergance of the echo 
canceller, it is desirable for the time delay to be determined sufficiently accurately 
with minimal computational complexity. 

Summary 

A time delay estimation system is provided for use with 
telecommunications systems to provide an estimate of signal delay. The time 
delay estimate system, according to a first embodiment, computes an average 
magnitude difference between energy envelope estimates of an original signal and 
a delayed signal. The time delay is determined according to a minimum in the 


average magnitude difference between the original and delayed signals. 

An absolute value generator receives an original signal and a delayed 
version of the original signal. The absolute value generator produces output 
signals representing the absolute value of the original signal and the delayed 
5 signal. A low pass filter is connected to the absolute value generator and receives 
the output signals from the absolute value generator. For each signal so received, 
the low pass filter substantially removes frequency components above a 
preselected frequency threshold to yield an estimate of the signal envelope. In that 
way, the low pass filter outputs an signal envelope estimate of the absolute value 

10 of the original signal and the delayed signal. To reduce computational 

requirements, a sample reduction device such as a decimation filter receives the 
filtered original and delayed signals to provide signals having a reduced sampling 
rate. A first buffer is connected to the sample reduction device and receives and 
stores the reduced original signal. A second buffer is connected to the sample 

15 reduction device and receives and stores the reduced delayed signal. Each buffer 
is connected to a difference processor which uses the information stored in the 
buffers to determine a series of differences between the buffered signals in order 
to determine an index producing a minimum difference. The time delay 
corresponding to the determined index then corresponds to the time delay 

20 between the original signal and the delayed signal. 

In a second embodiment, the difference processor is replace by a 
correlation processor. Both embodiments benefit from the robustness to noise 
design of the overall invention and the reduced computational complexity afforded 
by the data reduction. 
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Brief Description Of The Drawings 

The foregoing summary and the following detailed description of the 
preferred embodiments of the present invention will be best understood when read 
in conjunction with the appended drawings, in which: 

Figure 1 is a functional block diagram showing a first embodiment of the 
delay estimation system of the present invention; and 

Figure 2 is a functional block diagram showing a second embodiment of 
the delay estimation system of the present invention. 

Detailed Description 

Referring now to the drawings, wherein like numerals indicate like 
elements, Figs. 1 and 2 illustrate two embodiments of a delay estimation system 
10, 30 for estimating the delay introduced into a signal by a telecommunication 
system 12. The telecommunication system 12 receives an original signal OS and 
produces a delayed signal DS. Examples of telecommunications 12 may include 
any telecommunication system having an echo path causing such a delay between 
OS and DS. The delay estimation systems 10 includes an absolute value generator 
14, a low pass filter 16, a decimation filter 18, two buffers 20, 22, and a difference 
processor 24. 

The absolute value generator 14 receives the delayed signal DS and 
original signal OS, and produces an absolute value signal AVDS of the delayed 
signal and an absolute value signal AVOS of the original signal. These absolute 
value signals provide an estimate of the energy of each signal. A low pass filter 16 
is connected to the absolute value generator 14 for receiving the AVDS signal and 


the AVOS signal. The frequency response characteristic of the low pass filter 16 
is designed to substantially attenuate, or reject, a range of frequency components 
above a selected frequency. The low pass filter 16 alters the absolute value of the 
delayed AVDS to provide a filtered delayed signal FDS, which represents an 

5 estimate of the energy envelope of the AVDS signal. Likewise, the low pass filter 
is applied to the absolute value of the original signal AVOS and provides a filtered 
original signal FOS representative of the energy envelope of the AVOS signal. It 
is preferable to remove high frequency components to provide anti-aliasing for the 
decimation filter 18. Accordingly, the selected frequency is chosen with regard to 

10 the sampling rate of the sample reduction device 18. 

The decimation filter 18 is connected to the low pass filter 16 and receives 
the filtered delayed signal FDS and the filtered original signal FOS. The 
decimation filter 18 reduces the sampling rate of the filtered delayed signal FDS 
and the filtered original signal FOS, providing a reduced delayed signal RDS and 

15 reduced original signal ROS, respectively. The decimation filter 18 thereby 

reduces the computational requirements encountered in subsequent processing. 
The decimation filter 18 reduces the sampling rate by arbitrary factor M. 
Consequently, the decimation function reduces the subsequent computation 
requirements by a factor of M. The reduced delayed signal RDS is stored in a 

20 first buffer 20, and the reduced original signal ROS is stored in a second buffer 22, 
each buffer being connected to a respective output of the sample reduction device 
18. Each buffer stores a preselected number of samples, N. A continual estimate 
of the time delay can be obtained by removing the oldest sample from each buffer 
and adding a new sample to each buffer. Alternately, the entire buffer contents, 


may be periodically refreshed with new samples. 

A difference processor 24 is connected to receive the data stored in the 
first buffer 20 and the second buffer 22. The difference processor 24 performs an 
average magnitude difference calculation (AMDF) according to the following 
5 formula: 

Y N-l 

AMDF(i) = — 2k(*) " + 0| (1) 

where N is the number of samples stored in each buffer 20, 22, x^n) are the 
elements of the second buffer 22, x 2 (n) are the elements of the first buffer 20, and i 
represents a temporal shift between the data in two buffers 20, 22. For each value 
10 of i, the AMDF function returns a value representing the sum of the absolute value 
of the differences between the data of shifted first buffer 20, x 2 , and the second 
buffer 22, x v Optionally, the AMDF may be computed without dividing the sum 
by N. Thus, the AMDF reduces the computation requirements by eliminating 
multiplications. 

15 The value of the shift, i', yielding the smallest value of the AMDF 

represents an estimate of the delay between the original signal OS and the delayed 
signal DS. Specifically, the delay introduced by telecommunication system 12 
creates a commensurate offset between the data stored in the two buffers 20, 22. 
If the delay were to correspond to a precise integral value of i, i', the AMDF (i') 

20 would have a value of 0. However, the delay introduced to the delayed signal DS 
may not equate to a integral value of i, and the minimum value of the AMDF 
would be greater than 0. In this case, the value of i minimizing the AMDF would 
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provide an estimate of the delay between the original signal OS and the delayed 
signal DS. For larger values of the data reduction factor M, and thus fewer 
samples present in the buffers 20, 22, the larger the discrepancy is apt to be 
between the estimated delay time determined by the delay estimation device 24 

5 and the true delay time. Specifically, the time delay may be estimated by: 
td = (1/fs) M i', (2) 
where td is the delay time, fs is the sampling rate, M the data reduction factor, and 
i value that minimizes the AMDF. For example, for fs = 8000 samples per 
second, M = 8, if the minimal value of the AMDF occurs when i' = 4, the time 

10 delay is estimated to be 4 msec. This is termed the full bandwidth average 
magnitude difference. 

In another preferred embodiment of the AMDF, the computation 
requirements are further reduced by using a subset (window) of the reduced signal 
stored in the two buffers 20, 22. In this form the windowed AMDF is computed 

15 as 

1 B ~ l \ t 
AMDF (i) = — - x i( n + 0| (3) 

where B is less than N and represents the length of the window. Computation of 
the estimated delay for the windowed AMDF is performed in a manner analogous 
to the full bandwidth AMDF discussed above. 
20 In an alternative embodiment, the difference processor 24 may be 

replaced by a correlation processor to calculate a correlation function: 

^ N-l 

R(i) = Y^ in)x ^ n + i) (4) 

" n=0 
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where the variables have the same meaning as discussed above in connection with 
the AMDF. The time delay is estimated by a maximum value of the correlation 
function. Once the value of i' is determined that maximizes the value of the 
correlation function, computation of the time delay is given by equation (2) above. 
In addition, a windowing scheme may be applied to the correlation function in a 
similar manner to the windowed AMDF: 

«(i) = 7E x iW^(^o < 5 > 

where B is less than N and represents the length of the window. 

Further reduction in the computational requirements may be achieved by 
applying the decimation filter at an earlier stage in the time delay estimation 
process, as shown by the delay estimation device 30 in Fig. 2. The delay 
estimation system 30 includes a decimation filter 18' for receiving the delayed 
signal DS and original signal OS. The decimation filter 18' reduces the sampling 
rate of the delayed signal DS and the original signal OS, providing a reduced 
delayed signal RDS' and reduced original signal ROS', respectively. The 
decimation filter 18' thereby reduces the computational requirements encountered 
in subsequent processing. 

The reduced signals RDS 7 and ROS' are altered by an absolute value 
generator 14' to produce the absolute value of the reduced delayed signal AVDS' 
and an absolute value of the reduced original signal AVOS', respectively. A low 
pass filter 16' is connected to the absolute value generator 14' for receiving the 
AVDS' signal and the AVOS' signal. The frequency response characteristic of 


low pass filter 16' is designed to substantially attenuate, or reject, a range of 
frequency components above a selected frequency. The low pass filter 16' 
provides a filtered delayed signal FDS' and a filtered original signal FOS', which 
represent estimates of the envelopes of the AVDS' signal and AVOS' signal, 

5 respectively. Since sample reduction has already been preformed by the 
decimation filter 18', choice of the selected frequency may be made without 
concern of aliasing by the decimation filter 18'. 

The filtered delayed signal FDS' is stored in a first buffer 20, and the 
filtered original signal FOS' is stored in a second buffer 22, each buffer being 

10 connected to a respective output of the low pass filter 16'. Each buffer 20, 22 
stores a preselected number of samples, N. 

A difference processor 24 receives the data stored in the first buffer 20 
and the second buffer 22, The difference processor 24 may comprise an average 
magnitude difference function (AMDF) or may be replaced by a correlation 

15 processor as discussed above in connection with Fig. 1. In addition, the 

windowing techniques for the AMDF and correlation function may be applied to 
the delay estimation device 30. 

It will be recognized by those skilled in the art that changes or 
modifications may be made to the above-described embodiments without 

20 departing from the broad inventive concepts of the invention. For example, the 

original signal may be an impulse, so that the output of the difference processor 24 
provides an estimate of the delay time assisted with the impulse response function 
of the telecommunication system 12. It should therefore be understood that this 
invention is not limited to the particular embodiments described herein, but is 


Mlimmimmmmim 


intended to include all changes and modifications that are within the scope and 
spirit of the invention as set forth in the claims. 
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